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ABSTRACT 

We show that space-based telescopes such as the proposed Terrestrial Planet 
Finder Coronagraph will be able to detect the light scattered by the interstellar 
grains along lines of sight passing near stars in our Galaxy. The relative flux of the 
scattered light within one arcsecond of a star at 100 pc in a uniform interstellar 
medium of 0.1 H atoms cm~^ is about 10~^. The halo increases in strength with 
the distance to the star and is unlikely to limit the coronagraphic detection of 
planets around the nearest stars. Grains passing within 100 AU of Sun-like stars 
are deflected by radiation, gravity and magnetic forces, leading to features in the 
scattered light that can potentially reveal the strength of the stellar wind, the 
orientation of the stellar magnetic field and the relative motion between the star 
and the surrounding interstellar medium. 

Subject headings: ISM: dust — scattering — radiative transfer — techniques: 
high angular resolution 



1. INTRODUCTION 



Scattered starlight is a valuable tool for measuring the properties of interstellar dust. 
In opti cal reflection nebulae (IWitt et al.1 19921: ICalzetti et al.lll995l : iBurgh et al.ll2002r). dark 



clouds (IFitzgerald et a 



199ll : iMurthv fc Henrv 



II I 

.19761 : IWitt et a l. 1990) and the diffuse Galactic light (iHurwitz et al. 

I995I : IWitt et al.lll997l) the scattering grains typically lie many arc 



seconds from the illuminating sources. The dust can be nearer to or further from us than the 
source and the detected photons have been deflected through a wide range of angles. Fur- 
ther information about interstellar grain s comes from stellar X-rays scatt ered by foreground 



dust within a few arcseconds of the star (IMathis et al.lll995l : ISmith &: Dwek,1998( : IWitt et al. 
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200ll : iDraine fc Tarul2003l ). The measurements sample the whole column of material between 



us and the star because the grains are strongly forward scattering at X-ray energies. 

Here we show that high-contrast imaging with a space-based telescope such as the 
proposed Terrestrial Planet Finder Coronagraph enables a new kind of measurement for 
probing dust at projected separations similar to the X-ray measurements. The technique 
samples a range of scattering angles because at visible wavelengths the scattering is only 
mildly biased in the forward direction. As a consequence, most of the light comes from 
foreground grains lying close to the star. The scattered light halo carries information on the 
interstellar dust passing through the cavity opened in the interstellar gas by the stellar wind. 

Evidence for the existence of interstellar grains within our own heliosphere was obtained 
by the Ulysses spacecraft when 55 dust impacts as measured by the Cosmi c Dust Experimen t 



were identified as interstellar by their speed, mass and arrival direction (iGriin et al.lll994j ). 
The overall contribution of interstellar dust to the zodiacal cloud — the debris disk of the 
solar system — is unknown, although the fraction is presumably great er in the outer sola r 



system. Collisional debris from the asteroid belt dominates inside 3 AU (IGrogan et al.ll200ll ). 
However the interstellar grains, charged by photoionization and strongly influenced by the 
solar gravity, radiation pressure and magnetic field, are able to penetrate deep into the inner 
solar system. A uniform incoming spatial distribution of interstellar grains becomes strongly 
clumped as a function of particle size and phase of the solar cycle, although the contribution 
of this component to the all-sky thermal emission as vie wed from 1 AU is ne gligible due to 



the local dominance of asteroidal and cometary material (IGrogan et al.lll996l ) 



Starlight scattered by interstellar dust also is a source of noise in coronagraphic planet 
searches. We calculate the expected brightness and distribution of the scattered light and 
show that the halo is unlikely to affect the direct detection of Earth-like planets around the 
nearest stars. We compute the trajectories of interstellar grains under the stellar gravity, 
radiation and magnetic forces and show that the stellar wind produces observable signatures 
in the scattered light. For a large number of stars at distances greater than 100 pc, high- 
contrast imaging can potentially yield detailed information about the stellar wind and the 
adjacent interstellar medium. 

The remaining five sections of the paper cover the radiative transfer method (§2), the 
results for uniformly-distributed dust, obtained analytically assuming isotropic scattering 
(§3) and numerically including the anisotropy (§4), the results for dust passing through a 
model stellar wind (§5), and a summary and conclusions (§6). 
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2. RADIATIVE TRANSFER 

The transfer of visible light through the optically-thin interstellar medium is dominated 
by scattering and the radiative transfer equation at wavelength A is 

^ = ^ M^' ^ n)h{n')dn', (i) 

where Ix is the monochromatic specific intensity, tx the optical depth, ux the albedo of the 
scattering grains and $a the phase function describing scattering from all directions O' into 
the beam direction ft. The star is approximated by a disk of uniform specific intensity 
subtending a small angle 2il) as seen from a dust grain on our line of sight (figured]). The 
scattering phase function is taken to be constant across the face of the star, and the solid 
angle filled by the star is tt?/'^. The transfer equation reduces to 

^ = -cox<^xi^)7i^'ll (2) 

dTx 

where (j) is the scattering angle between the rays joining star to grain and grain to observer. 
The intensity of the scattered light is calculated by integrating eq. [2] along the hne of sight. 
The angles and ip vary with the distance from the observer, z. Writing the distance in place 
of the optical depth using drx = Xxpdz, neglecting emission from any distant background, 
and taking the simple case of an interstellar medium with uniform opacity, the intensity is 

POO 

h = XxpuJxIl / $A(0)vrV'rf^. (3) 

^0 



3. ISOTROPIC SCATTERING 

We first solve equation [3] assuming isotropic scattering, with $a = 1/ (4vr) independent 
of (f) and z. Along a line of sight offset by a small angle 9 from the direction to the star, the 
stellar angular radius ip varies as 

^{^) = (4) 

where the star has radius and lies at distance D. With a change of variable to y = 
(1 — z/D)/9, the solution to the transfer equation is 

Ix{e) xxP^xrl f'/' dy 

(5) 



/* ADe 1 + 1/2 



observer 



D 



star 



Fig. 1. — Sketch of the scattering geometry for a star at distance D. A representative interstellar 
grain lies along a line of sight offset from the star by an angle 6 and is at a distance z. The star 
subtends an angle 2il) as seen from the grain. Stellar photons reach our line of sight on scattering 
through an angle </>. 



The integral is very nearly vr because 1/^ is almost 
with y. 



-oo and the integrand falls off rapidly 



As a specific case we take the Sun observed from a distance of 100 pc at wavelength 
547 nm through a uniform interstell ar medium with 0.1 H atoms cm~^, similar to the de nsities 
found in the solar neighborhood (iBruhweiler fc Kondd Il982l : ISlavin fc FrischI l2002l ). For 
Milky Way grains, the total opacity y\p = 4.9 x 10~^^ cm^ H"^ x 0.1 H cm"^ = 4.9 x 10~^^ 
cm~^ and albedo uj\ = 0.66 (lDraindl2003l ). The optical depth to the star X\pD ~ 10^^. The 
predicted surface brightness ratio at an angular separation 6 is 



h{0) TTXxpuJxrl 



we 



8x10-^^ 

VO.l cm 



100pc\ /O.l" 



and falls off inversely with the separation. The integrated flux within an angle 9^ 
proportional to the optical depth, 




(6) 



IS 



(7) 



and the flux within one arcsecond is eighty times greater than the nominal threshold 10 F; 
for detection by the proposed Terrestrial Planet Finder Coronagraph. 
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4. ANISOTROPIC SCATTERING 

Interstellar grains scatter visible light anisotropically. We treat this effect using a stan- 
dard model for grains in the diffuse interstell ar medium of the Milky Way, with the phase 



function $a(0) computed from Mie theory by iDraind (120031 ) . Forward scattering is moder- 



ately preferred and at the chosen wavelength of 547 nm, each event has a 50% chance of 
deflecting the photon through an angle less than 43°. The interstellar medium is uniform 
with the same density as in section [31 The transfer equation [2] is numerically integrated 
along the line of sight to the observer using a fourth-order Runge-Kutta method. We have 
checked that the method reproduces the analytic solution given by equation in the case of 
isotropic scattering. Including the preferential forward scattering yields a halo 48% brighter. 
The results for a solar analog star at a distance of 100 pc are shown in figure El Also shown 



for comparison is a solar analog exozodiacal cloud from appendix l.B of iLevine et al.l (120061 ). 
The light scattered from the interstellar grains is brighter than that from the exozodiacal 
cloud outside 50 milliarcsec. The exozodiacal cloud is inclined 60° from face-on and its 
brightness was calculated including anisotropic scattering. 

Much of the interstellar scattered light comes from grains passing close by the star, 
due to the inverse square falloff of the stellar illumination. Most comes from grains whose 
distance from the star is less than ten times the projected separation in the plane of the sky. 
The preferential forward scattering has a smaller effect, making the dust in front of the star 
contribute more than the dust at the same distance behind. The contributions from dust 
at different positions along a line of sight passing 0.1 arcsecond from a star at 100 pc are 
shown in figure [3l Overall, the scattered light halo is most sensitive to the distribution of 
interstellar dust in the vicinity of the illuminating star. 

Measurements and models of the optical properties of the dust in the diffuse int erstellar 
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the albedo and the scattering asymmetry as measured by the mean cosine scattering angle 
g = (cos 0) at 547 nm. The appearance of the halo is qualitatively independent of the dust 
properties over the allowed range. The halo brightness is simply proportional to the albedo 
(equation [6|). The scattering asymmetry affects both the halo brightness and the location 
of the grains giving most of the light. We examine the dependence by re-computing the 
halo using a standard Henyey-Greenstein phase function with asymmetry parameters up to 
the maximum g = 0.8 consistent with the range of the measurements. Half the photons are 
defl ected through angles less than 19°. The resulting halo is 2.19 times brighter than with 



the iDraind (120031 ) phase function. The fraction of the light 0.1 arcsecond from the star that 



is scattered by grains lying more than 200 AU from the star rises from 7% to 21%. A star 



- 6 - 




Fig. 2. — Surface brightness versus separation from the star for a solar analog exozodiacal cloud 
(dashed curve) and uniformly-distribut ed interstellar g rains, scattering either isotropically (dotted 
line) or anisotropically according to the lDraind (120031 ) phase function (solid line). A crossed circle 
marks the approximate angular separation and surface brightness of an exo-Earth at t he nominal 
angular resolution of the Terrestrial Planet Finder Coronagraph ( Levine et al.ll2006l ). The star 
lies at a distance of 100 pc. 



surrounded by the more forward-throwing grains has a brighter halo with a shghtly lower 
contrast for any features arising from dust passing close by the star. 



PROBING ASTROSPHERES 



Sub-micron interstellar grains entering the solar syster n are deflected by the solar g ravity, 
radiation pressure and solar wind electromagnetic forces (I Gust af son &: Misconilll979l ). The 
highest densities of grains 0.1 micron in radius occur where the particles pile up against the 
solar wind and where the magnetic f orces concentrate particles over the solar magnetic poles 
(iGrogan et al.lll996l : lLandgrafll2000l ). In this section we explore whether similar effects can 
be detected in the scattered light halos around other stars. 
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Fig. 3. — Contribution to the surface brightness from grains at different locations along the line 
of sight passing 100 mas or 10 AU from the star in figure [21 The three panels show isotropic 
scattering (top) and anisotropic scattering by grains in front of (center) and behind the plane 
containing the star (bottom). The horizontal axes give the distance from the star plane in AU 
(above the chart) and parsecs (below). Most of the light is scattered into our line of sight by 
grains within 100 AU of the star. 



5.1. Stellar Wind Model 



Th e star is taken t o be a solar analog and its w ind is modeled foUowing lGustafson fc Misconi 



fll979h . iGrogan et all fll996f ) and 



andgraj (l2000f) by applying the expanding corona model 



of the solar wind magnetic field (iParkerl 119581 ). The gas streaming outward from the star 



draws out the stellar magnetic field lines so that near the star the field is approximately 
radial. The more distant field lines are dragged into the shape of an Archimedean spiral as 
the star rotates. The radial and azimuthal components of the field are 

Br = Bo/r^ (8) 



and 



Ba 



-BoUjb'^{r — h) sin a/Kj^r^ 



(9) 



where u is the stellar angular velocity, a is the astrocentric colatitude angle, Vsw is the stellar 
wind speed (chosen as 400 km s~^) and h is the distance at which the reference field strength 
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Bq is taken. The Lorentz force is then 



[{Vg xB)- {Vs^ X B)] 
where q is the grain charge and Vg the grain speed. The grain charge q is given by 



(10) 



q = AneoUs (11) 

where s is the grain radius, Eq = 8.86 x 10~^^ CV~^m~^ is the permittivity of free space and 
U is the grain surface potential. Dust in the stellar wind cavity is expected to be positively 
charged due to the dominance of stell ar photoioniz ation over the competing mechanism of 
the 'sticking' of stellar wind electrons (jMukailll98ll ). Since the stellar wind density roughly 
follows an inverse square law decrease with astrocentric distance, particles are charged to 
+5-10 volts with little distance variation. We adopt a surface potential U of +5 volts. 
Notice that the dominant electric part in equation [10] is proportional to B^ and therefore 
to r~^. This can make the instantaneous Lorentz force domin a te on small grains at large 
astrocentric distances. The ballerina model proposed by lAlfvenI (119771 ) gives an explanation 
for solar system spacecraft measurements which indicate that near the ecliptic plane the 
interplanetary magnetic field is directed inward in certain regions and outward in others. 
In this picture a current sheet separates plasma from either hemisphere carrying fields of 
opposite polarity. As the sheet rotates with the Sun the small up and down displacements 
in the sheet, similar to the wave motion of the skirt of a spinning ballerina, explain the 
observed effect. We use a simplified model of this phenomenon to describe the magnetic 
field strength averaged over a stellar rotation by the inclusion of a current sheet which 
changes its tilt angle e to the stellar equatorial plane at a con stant angular ve locity of 
(l/ll)7r radians per year over the stellar magnetic activity cycle (iGustafsonl Il994j ) . so that 
the sheet completes one 360° rotation in 22 years. The averaged field at latitude A is then 
a fraction 2 arcsin(tan74/ tane)/7r of the unipolar field. The line of nodes defined by the 
intersection of the current sheet and stellar equatorial plane will rotate through 360° in one 
27 day stellar rotation. The effect of radiation pressure is to reduce the stellar gravitational 
force, the magnitude of the radiation pressure being given by 

S 



pr 



-Q 



pr 



(12) 



where S is the flux, c is the speed of light and Qpr is the radiation pressure efficiency of the 
particle. 



5.2. Particle Properties 



Grains larger than one micron are rare in the diffuse interstellar medium according 
to extinction curve fitting results (IMathis. Rumpl &: Nordsieckl 119771 : iDraine fc Led Il984j : 
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Weingartner fc Draind l200ll ). Grains smaller than 0.01 micron have little optical cross- 
section and furthermore are largely excluded from the heliosphere by magnetic forces due to 
their high charge-to- mass ratios. Grains in the 0.1 micron size range experience a radiation 
pressure force larger t han the gravitation al force but nevertheless are able to penetrate deep 
into the solar system ( Grogan et al.lE996l ). Most of the scattered light from interstellar grains 
passing through the astrosphere of a solar analog star is likely due to grains about 0.1 micron 
in size. We therefore concentrate on grains with radii between 0.02 and 0.3 microns. 

The ratio (3 of the radiation pressure force to the gravitational force may be expressed 
in terms of Qpr as 



5.7 X lo-^g 



pr 



sp 



(13) 



where p is the density of the particle (iBurns et al.lll979l ). The optical efficiencies of the 



grains are found from Mie theo ry. We assume astronomical silicate composition and optical 
constants (iLaor fc Drainelll993l ). Clearly, the particle spatial distributions will vary with the 



particle size; for example, the gravitational force will dominate for larger particles, whereas 
the interaction with the magnetic field will become increasingly important at the small 
particle end of the size range. The distribution of particles will also depend on the stellar 
magnetic cycle, so that the particle positions will represent a 'snapshot' in time rather than 
an equilibrium configuration. 



5.3. Dynamical Calculations 

Taking into account gravity, radiation pressure and the stellar magnetic field model 
outlined above, we follow the trajectories of a uniform distribution of dust particles at 
infinity as they approach the star from the upstream direction at a speed of 26 km s~^. 
The equations of motion are integrated at small intervals of time and the position, speed 
and acceleration of the particles are calculated. The local properties of the stellar wind and 
magnetic field are updated as the integration progresses. We continue the integration until 
the heliosphere analog is filled with dust particles and stop the run at a given phase of the 
stellar magnetic field cycle with the magnetic field orientation matching that of the Sun in 
the year 1990. 

The dust density is then calculated by binning the particles in cubes 2 AU on a side 
and normalizing to the upstream interstellar density. The resulting distributions are broadly 
similar for particles 0.02 to 0.3 pm in radius. The distribution of 0.1 pm grains is shown in 
figure m The particles are pushed by the magnetic forces into the equatorial current sheet, 
and pile up 10 AU upstream from the star (left panel), where the density reaches six times 
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Fig. 4. — Interstellar dust density near a Sun-like star. The star is centered in each image and 
the dust enters from the right, as shown by the arrow. At left is a slice through the equatorial 
plane and at right is the polar plane containing the velocity vector. The images are 400 AU on 
a side. The color scale is linear between zero (black) and the maximum density (red). Focusing 
by the stellar gravity, radiation and magnetic forces enhances the density of the 0.1 fim particles 
up to 6 times at left and 15 times on the right. 

that in the interstellar medium. The greatest density enhancement of all is a factor fifteen 
and occurs in two fins over the magnetic poles and 60 AU downstream from the star (right 
panel). These result from the focusing of the grains by the inclined magnetic fields off the 
equatorial plane. 



5.4. Synthetic Images 



Synthetic images of the scattered starlight are constructed by solving the transfer equa- 
tion [2] to find the intensity on a grid of sky points, treating the anisotropic scattering as 
described in section H] with the iDraind (120031 ) Mie theory phase function. No exozodiacal 
cloud is included. Selected results are shown in figures [5] and [61 In addition to an overall 
deficit of scattered light relative to the case of uniformly-distributed dust, several features 
produced by the magnetic field are clearly visible, including the density enhancements in the 
equatorial disk and the polar fins. 



6. CONCLUSIONS 



We have shown that stars appear surrounded by halos of light scattered by interstellar 
dust. The halos are bright enough for detection with the proposed Terrestrial Planet Finder 
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Fig. 5. — Model images of the starlight scattered by 0.1 fim interstellar grains around a Sun-like 
star at a distance of 100 parsec. The star is viewed looking upstream (left panel) and from the 
side (right panel). The cross and arrow indicate the flow of the approaching grains. The images 
are 2 arcsec or 200 AU across, and were divided by the profile for uniformly-distributed dust. Dark 
regions show a deficit of scattered light where magnetic and radiation forces exclude particles 
near the star. The gray scale is linear between the minimum brightness ratio of about 0.1 (black) 
and the maximum of unity (white). 



Coronagraph around stars at distances of 100 pc and greater. The halo reveals the distri- 
bution of dust in the vicinity of the star, providing opportunities to probe the interstellar 
medium and measure the properties of the stellar wind. Particles of differ ent sizes hav e 



differ ent distributions, with larger particles approaching the star most closely (IGrogan et al. 



19961 ). This size-sorting can potentially be used to measure the distribution of grain sizes at 



selected locations in the interstellar medium and to resolve whether particles larger than one 
micron are present in significant nu mbers, as suggested by dust impact detector da ta from 
the Ulysses and Galileo spacecraft (IFrisch et al.lll999l : IWeingartner fc Drainell200ll ). Since 
the grains are deflected by stellar wind magnetic forces, the features in the scattered light 
halo can be used to estimate the wind speed and the relative motion between the star and 
the surrounding interstellar medium. The results will complement the wind mass loss rates 
estimated using th e Lyman-a absorption by the gas collected near stellar wind bow shocks 
(jWood et al.l 120051 ). The scattering halos will make it possible in addition to determine 
the wind magnetic field strength and orientation (figure [5]) and, by obse rvations spanning a 
decade or so, the variation over the stellar activity cycle (jLandgraill2000l ) . Imaging the halos 
will yield fresh information about the winds of sta rs with a range of ages and could lead 
to a better understanding of stellar spin histories (IBouvier et al.l 119971 ) and the long-term 
evolution of stellar activity. 



This work was carried out at the Jet Propulsion Laboratory, California Institute of 
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Fig. 6. — Model profiles of the starlight scattered by 0.1 /im interstellar grains around a Sun-like 
star at a distance of 100 parsec. The solid line is for uniformly-distributed dust. The dashed 
line shows the effects of stellar gravity, radiation and magnetic forces pushing aside dust near 
the star, and is a horizontal cut through the center of figure O left panel. The vertical axis is 
the count rate per pixel on a logarithmic scale, in units of the stellar count rate. The pixels are 
20 milliarcsec on a side. 

Technology with support from the internal Research and Technology Development Fund 
Innovative Spontaneous Concepts Program. 
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